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ABSTRACT 

We use a homogeneous sample of gamma-ray bursts (GRB) extracted from 5.9 years 
of BATSE DISCLA data (Schmidt 1999) and a variety of assumed broken power-law lu- 
minosity functions to derive GRB luminosities and space densities. Luminosity functions 
that are narrow or exhibit no density evolution produce expected redshift distributions 
that are incompatible with the observation of a GRB redshift of 3.4. For qo = Q.l and 
density evolution rising to 10 at z = 1, we find for a variety of slopes of the luminosity 
function values of the local space density around 0.18 Gpc~^ y~^ with a range of only 
40%. Characteristic 50 — 300 keV peak luminosities exhibit a range of a factor of 6 
around 6 x 10^^ erg s~^, corresponding to a characteristic total luminosity of 1.2 x 10^^ 
erg in the 10 — 1000 keV band. For Qq = 0.5, densities are higher and luminosities lower, 
both by a factor of 2.5. The local emissivity of GRBs is 1.0 x lO^'^ erg Gpc^^ y~^ in 
the 10 - 1000 keV band. 

Subject headings: cosmology: observations — gamma rays: bursts 

1. Introduction 

Soon after the launch of the Compton Gamma Ray Observatory in 1991, observations of 
gamma-ray bursts (GRB) with the Burst and Transient Source Experiment (BATSE) established 
two important global properties of GRBs: an isotropic sky distribution (Meegan et al. 1992a) and a 
radial distribution incompatible with a uniform population in euclidean space (Meegan et al. 1992b). 
The evidence for a non- uniform radial distribution was based on either the N(> P) distribution or 
the V/Vmax distribution of the observed GRBs. The N{> P) distribution of the peak flux P should 
be a power law of slope —3/2 for a uniform population in euclidean space. The observed distribu- 
tion was shallower, particularly at lower fluxes (Meegan et al. 1992a). For a uniform population, 
the V/Vmax distribution would be uniform between and 1, with < V/Vmax >= 0.5 (Schmidt, 
Higdon and Hueter 1988). Early BATSE observations yielded < V/Vmax >= 0.35, significantly 
different from the value expected for homogeneity (Meegan et al. 1992b). 

A cosmological distance scale of GRBs is compatible with the isotropic sky distribution and 
the radial distribution (Paczynski 1992). This has been confirmed by the first observation of a large 
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redshift for a GRB afterglow (Metzger et al. 1997). The distance scale in this cosmological scenario 
for GRBs has been derived from the observed A^(> P) distribution, usually on the assumption that 
GRBs are standard candles (Pendleton et al. 1996; Wijers et al. 1998; Totani 1999). Krumholz, 
Thorsett, and Harrison (1998) have shown that relaxing the standard candle assumption renders a 
broad range of models consistent with the BATSE N{> P) relation. 

In this Letter, we are reporting the results of setting the cosmological distance scale of GRBs, 
assuming a full luminosity function and using the euclidean < V/Vmax > value as a distance 
indicator. We use the BD2 sample of GRBs derived from 5.9 years of BATSE DISCLA data, 
briefly described in Sec. 2. The characterization of the luminosity function and its evolution, and 
the derivation of predicted distributions of luminosity, V/Vmax, redshift, and flux are discussed in 
Sec. 3. The results for a variety of luminosity functions are presented in Sec. 4 and the conclusions 
are summarized in Sec. 5. We are assuming in this paper a Hubble constant Hg = 70 km s~^ 
Mpc~^ and zero cosmological constant. 

2. The BD2 Sample of Gamma-Ray Bursts 

We use a large homogeneous sample (the BD2 sample) of 1391 GRBs derived from BATSE 
DISCLA data, consisting of the continuous data stream from the eight BATSE detectors in four 
energy channels on a time scale of 1024 msec (Fishman et al. 1989). This sample is a revision, 
described below, of the BDl sample which resulted from a search of DISCLA data over the period 
TJD 8365 — 10528 (Schmidt 1999). The BDl sample was based on a trigger algorithm that used 
the background both before and after the onset of the burst and required an excess of at least 5a 
over background in at least two detectors in the energy range 50 — 300 keV. 

In the process of classifying triggers for the creation of the BDl sample (Schmidt 1999), we had 
accepted 1018 DISCLA triggers that were within 230 sec of a GRB in the BATSE catalog (Meegan 
et al. 1999) as genuine GRBs. In addition, we classified another 404 DISCLA triggers as GRBs. 
In a subsequent revision, we have now inspected the output of the BATSE detectors over a time 
interval of 12,000 sec around each of the 404 GRBs not in the BATSE catalog, as well as those in 
the catalog for which times or positions differed appreciably. In the process, we rejected 7 triggers 
as caused by source fluctuations, 18 turned out to be parts of other GRBs of long duration, and 
6 were identified as the soft repeater SGR1806-20 and rejected. As a consequence of the revision, 
the BD2 sample now contains 1391 GRBs, of which 1013 are listed in the BATSE catalog, another 
377 are classified as GRBs, and one as a probable GRB. 

For the purpose of this paper, we can characterize the BD2 sample as follows. The number of 
GRBs is 1391. The sample effectively represents 2.003 years of full sky coverage (Schmidt 1999), 
so the rate is 694 GRBs per year. The average euclidean V/Vmax is 0.334 it 0.008. The limiting 
flux has a distribution G{Piim) that has been derived as follows. For 104 positions of fixed celestial 
coordinates distributed isotropically around the sky, we checked every 100 sec during every tenth 
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day whether the sky position was above the horizon and whether the time did not fall in an exclusion 

window set up to avoid interference or bad data (Schmidt 1999). If so, we converted the limiting 
count of 5 times the square root of the observed background in the second brightest illuminated 
detector into a limiting flux Pum- The resulting distribution G{Piim) has a median at 0.37 ph cm~^ 
s~^, and 10 and 90 percentiles at 0.29 ph cm~^ s~^ and 0.51 ph cm~^ s"^, respectively, in the 
50 - 300 keV band. 



3. Derivation of the Distance Scale 

Since the number of redshifts of GRBs is as yet too small for a statistical derivation of the 
luminosity function, we will instead assume a full luminosity function and then use the observed 
number of GRBs and their euclidean < V/Vmax > in the BD2 sample to derive properties such 
as the local space density and the characteristic luminosity L*. We will find that these properties 
vary relatively little as we change the shape of the luminosity function. 



3.1. The Luminosity Function 

Experience has shown that the differential luminosity function of many types of extragalactic 
objects can be represented as a broken power law (cf. e.g. Hasinger (1998) for X-ray active 
galactic nuclei). We will use a broken power law, and in addition introduce upper and lower limits 
of luminosity, thus allowing both narrow (standard candle) and more realistic broad luminosity 
functions. We will also assume density evolution, to be introduced below. 

The local luminosity function of peak ORB luminosities L, defined as the co-moving space 
density of GRBs in the interval log L to log L + d log L, is 

$o(L) = 0, for logL < logL* - Ai, (la) 

$„(L) = Co{L/L*r\ for logL*-Ai <logL<logL*, (lb) 

$„(L) = Co{L/L*r\ for logL* <logL<logL* + A2, (Ic) 

$o(L) = 0, for logL > logL* + A2. (Id) 



Given that most types of extragalactic objects show evolution, and considering that the rate 
of star formation may be relevant for GRBs (Madau, Pozetti, and Dickinson 1998), we introduce 
density evolution p{z). In most cases, we have assumed 

p{z) = {l + zf, for 0<z<zp, (2) 
p{z) = {l + ZpY, for z>zp. (3) 

While any type of evolution can be explored (Che, Yang, and Nemiroff 1999), we have used in this 



-4- 



Letter only the case Zp = 1.0 and p{1.0) = 10 without a dechne toward larger redshifts (Steidel et 
al. 1999). The luminosity function at redshift z is ^ziL) = ^o{L')p{z). 



3.2. Predicting the Sample Properties from the Luminosity Function 

In principle, with the full characterization of the BD2 sample given in Sec. 2 and the luminosity 
function ^z{L) with all its free parameters set, we can predict or model all properties of the sample. 
Our procedure will be to assume values for the parameters Ai, A2, ai, and 02, and then to vary 
L* until the predicted value of the euclidean < V/Vmax > equals the observed value in the BD2 
sample, and to set Cg to fit the observed number of GRBs. 

The modeling procedure involves the derivation of the peak flux P{L, z) of a GRB of peak 
luminosity L observed at redshift z (Schmidt and Green 1986), 

where A{z) is the bolometric luminosity distance and C{Ei,E2) is the integral of the spectral 
energy distribution between Ei = 50 keV and E2 = 300 keV. We use a Band spectrum (Band et 
al. 1993) for the energy distribution with a = —1.0, /? = —2.0, and Eg = 200 keV in the restframe. 

Objects with luminosity L observed in a part of the BD2 sample with flux limit Pn^ are 
detectable to a maximum redshift Zjnax{L,Piim) that is easily derived from eq.(4). The total 
number of objects in the sample is 

/r rZmaxiLjPlim.) 
$o(L) d log L y G(Pam) dPii^ p{z) (dV/dz) dz (5) 

Based on this model we can derive the distributions of the peak flux P, the euclidean V/Vmax = 
{P/Piim)~^^'^, the peak luminosity L, and the redshift z. These form the basis for the discussion 
of results below. 



4. Results for a Variety of Luminosity Functions 

Wc exhibit in Table 1 numerical results for 18 luminosity function models, nine each for the 
cases Qo = 0.1 and Qo = 0.5. All are based on a Hubble constant Ho = 70 km s^^ Mpc~^ and zero 
cosmological constant. Among the nine cases, two have zero density evolution, six have evolution 
as described in eqs. (2) and (3), and one is an exponential of cosmic time. In the six cases, we 
explored variations of the extent of the luminosity function below L*, as well as of the logarithmic 
slopes of both parts, cf. Table 1. 

As an example, model 14 uses the deceleration parameter Qo = 0.1 and logarithmic slopes of 
—0.5 and —2.0 for the power law luminosity function, which extends from L*/10 to 100 L*. The 
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model has density evolution p{z) = (1 + z)^'^^ for < z < 1 and p{z) = 10.0 for z > I. Figure 1 
shows the relation between the model values of < V/Vmax > and logL*. The observed value of 
< V/Vmax >= 0.334 lb 0.008 corresponds to logL* = 51.76 ib 0.08, where the error is the formal 
error corresponding to that in < V/Vmax >■ The value of Co is set by the total number of GRBs in 
the BD2 sample. With the luminosity function of model 14 fully characterized, we can derive the 
predicted distributions in the BD2 sample of the luminosities, the fluxes, the V/Vmax values and 
the redshifts; the last three are illustrated in Figure 2 together with the luminosity ftinction. Apart 
from the obvious onset of incompleteness in the observed fluxes, the predicted A^(> P) agrees well 
with the observed distribution. The model distribution of V/Vmax shows fairly good agreement 
with the observations. The distribution of redshifts is broad, with some redshifts expected as large 
as 6. The median redshift is 1.5. 

Besides the values of log L* and Cq derived for each model, we show in Table 1 also the local 
{z = 0) space density po of GRBs integrated over the full range of luminosities L, and the probability 
P3.4 that a GRB in the BD2 sample has a redshift z > 3.4. This probability is of interest, since one 
of the first observed GRB redshifts was 3.4 (Kulkarni et al. 1998). For the models with a standard 
candle luminosity function or with zero evolution (nos. 11 — 13 and 51 — 53), the values of P3.4 
are so low that they are inconsistent with the observed redshift of 3.4. Model 53 is similar to the 
standard candle luminosity models considered by Wijers et al. (1998) and Totani (1999). Wijers et 
al. (1998) derived a luminosity about 3.5 times larger, and a local density 2.5 times smaller than 
those given in Table 1, while Totani (1999) found a luminosity similar to our value. The broad 
luminosity functions considered by Krumholz, Thorsctt, and Harrison (1998) cover a range from 
10™ — 10^^ erg s~^, overlapping with those in models 54 — 58. 

We show in Figure 3 a plot of po vs. L* for all models. There are clear systematics that can 
be summarized as follows. From Qq = 0.1 to Qq = 0.5 local densities increase and characteristic 
luminosities decrease, both by a factor of about 2.5. The effect of the density evolution adopted 
is to decrease local densities by a factor of 15 — 20 and to increase L* by a factor of 3. For given 
cosmology and evolution, local densities range over a factor less than 2 and L* over a factor of 6. 
Given that po and L* tend to be inversely correlated, we derived for each model Eoy,t, the total 
peak output of GRBs per unit volume at z = 0. As shown in Table 1, Eo^t is remarkably stable 
at 5 X 10^° erg s~^Gpc~^ y~^ for the adopted evolution. For zero evolution, Eout is some 5 times 
larger. 

If the gamma radiation of GRBs is beamed or collimated, say over iirF steradians, then all 
luminosities should be multiplied by F, and all space densities divided by F. The total peak output 
Eout is unaffected by beaming. 
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5. Summary 

We have shown that using primarily < V/Vmax wc can derive properties of GRBs of useful 
accuracy. We have used redshift information only in arguing that the observation of a redshift of 
3.4 for a GRB is unlikely if the luminosity function is a standard candle, or if the GRB rate does 
not evolve with redshift. For density cvohition rising to a factor of 10 at z = 1 and Qq = 0.1 (for 
Qo = 0.5, cf. Table 1), the local space density is 0.18 Gpc^^ y~^, with a range of only 40% among 
the models considered. Characteristic peak luminosities range over a factor of 6 around 6 x 10^^ 
erg s^^ in the 50 — 300 keV band. 

Luminositics"'^ derived in this study are peak luminosities (strictly per 1024 msec), since the 
detection of GRBs in the BD2 sample is on that time scale. We estimate the total energy radiated 
by GRBs in gamma rays by integrating over the time profile and by extending the spectral range 
to 10 — 1000 keV. To account for the total energy radiated in the 50 — 300 keV band integrated 
over the duration of the burst, we consulted data in the BATSE catalog (Meegan et al. 1999). We 
find that for the 100 — 200 strongest bursts, the ratio of fluence over peak flux is around 9 — 10 s. 
The ratio of the energy radiated in a band of 10 — 1000 keV over that in the 50 — 300 keV range 
based on the Band spectrum (Band et al. 1993) with a = —1.0, /? = —2.0, and Eq = 200 keV is 
2.1. We conclude that total luminosities over the 10 — 1000 keV band are 20 times larger than the 
50 — 300 keV peak luminosities. Employing this factor, we obtain for Qo = 0.1 a characteristic GRB 
luminosity of 1.2 x 10^^ erg in the 10— 1000 kcV band. Similarly, the local cmissivity of GRBs Eo^t 
discussed in Sec. 4 is 1.0 x 10^^ erg Gpc~^ y~^ in the 10 — 1000 keV band with very little range 
among the models. 
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Fig. 1. — Plot of the euclidean value of < V/Vmax > versus peak luminosity for luminosity function 
model 14, cf. Table 1. 

Fig. 2. — Properties of luminosity function model 14, cf. Table 1. Panel (a) shows the luminosity 
function (Gpc~^ y"-*^ per unit log-L). The break luminosity L* has been set by V/Vmax and 

the normalization by the observed rate of GRBs. Panel (b) compares the model euclidean V/Vmax 
distribution for the BD2 sample with the observed distribution. Panel (c) compares the model 
flux distribution A'^(> P) for the BD2 sample with the observations. Panel (d) shows the expected 
distribution of redshifts for the BD2 sample, with an indication of the probability P3.4 of a redshift 
exceeding 3.4. 

Fig. 3. — Plot of the local GRB density versus peak luminosity L* for 18 luminosity function 
models, cf. Table 1. The systematic effects of the cosmological model {qo = 0.1 and Qo = 0.5) 
and of evolution are easily discerned. For a given cosmological model and evolution, local densities 
range over a factor of 2, and luminosities over a factor of 6 for the different luminosity function 
models considered. 
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Table 1. Properties of various luminosity functions models. 
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^In all cases we use Hq = 70 km s Mpc and zero cosmological constant. 

^L* is the peak luminosity in erg s~^ in the 50 — 300 keV band, if the GRB is radiating 
isotropically. 

■^Units are Gpc~^ y"^ per unit logL. 

'^Po is the local {z = 0) GRB rate, in units of Gpc~^ y~^. 

^Probability of z > 3.4 for a GRB in the BD2 sample. 

^Local peak energy output of GRBs in the 50 — 300 keV band, in units of 10^° erg s~^ 
Gpc-3 y-i. 

^In this case p{z) = e^"^^^^ where t{z) is the light travel time, expressed in the age of the 
universe. 
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